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Abstract 

Understanding the nature of dense particle packings is a subject of intense research in the 
physical, mathematical and biological sciences. The preponderance of previous work has focused 
on spherical particles, and very little is known about dense polyhedral packings. We formulate the 
problem of generating dense packings of nonover lapping, non-tiling polyhedra within an adaptive 
fundamental cell subject to periodic boundary conditions as an optimization problem, which we call 
the Adaptive Shrinking Cell (ASC) scheme. This novel optimization problem is solved here (using 
a variety of multi-particle initial configurations) to find the dense packings of each of the Platonic 
solids in three-dimensional Euclidean space M"^, except for the cube, which is the only Platonic solid 
that tiles space. We find the densest known packings of tetrahedra, icosahedra, dodecahedra, and 
octahedra with densities 0.823 . . ., 0.836 . . ., 0.904 . . ., and 0.947 . . ., respectively. It is noteworthy 
that the densest tetrahedral packing possesses no long-range order. Unlike the densest tetrahedral 
packing, which must not be a Bravais lattice packing, the densest packings of the other non-tiling 
Platonic solids that we obtain are their previously known optimal (Bravais) lattice packings. We 
also derive a simple upper bound on the maximal density of packings of congruent nonspherical 
particles, and apply it to Platonic solids, Archimedean solids, superballs and ellipsoids. Provided 
that what we term the "asphericity" (ratio of the circumradius to inradius) is sufficiently small, 
the upper bounds are relatively tight and thus close to the corresponding densities of the optimal 
lattice packings of the centrally symmetric Platonic and Archimedean solids. Our simulation 
results, rigorous upper bounds, and other theoretical arguments lead us to the conjecture that 
the densest packings of Platonic and Archimedean solids with central symmetry are given by 
their corresponding densest lattice packings. This can be regarded to be the analog of Kepler's 
sphere conjecture for these solids. The truncated tetrahedron is the only non-centrally symmetric 
Archimedean solid, the densest known packing of which is a non-lattice packing with density at least 
as high as 23/24 = 0.958333 .... We discuss the validity of our conjecture to packings of superballs, 
prisms and anti-prisms as well as to high-dimensional analogs of the Platonic solids. In addition, we 
conjecture that the optimal packing of any convex, congruent polyhedron without central symmetry 
generally is not a lattice packing. Finally, we discuss the possible applications and generalizations 
of the ASC scheme in the predicting the crystal structures of polyhedral nanoparticles and the 
study of random packings of hard polyhedra. 
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I. INTRODUCTION 



Particle packing problems are ancient, dating back to the dawn of civilization. Bernal 
has remarked that "heaps" (particle packings) were the first things that were ever measured 
in the form of basketfuls of grain for the purpose of trading or of collection of taxes 
Dense packings of hard particles have served as useful models to understand the structure 
of liquid, glassy and crystal states of matter 

aaa, 

granular media 



materials 



5j, and heterogeneous 



]. Understanding the symmetries and other mathematical properties of the 
densest packings in arbitrary dimensions is a problem of long-standing interest in discrete 
geometry and number theory [3]. 

A large collection of nonoverlapping solid objects (particles) in d-dimensional Euclidean 
space M"' is called a packing. The packing density cf) is defined as the fraction of space 
covered by the particles. A problem that has been a source of fascination to mathe- 
maticians and scientists for centuries is the determination of the densest arrangement (s) 
of particles that do not tile space and the associated maximal density (pmax Q]- Finding 
the maximal-density packing arrangements is directly relevant to understanding the struc- 
ture and properties of crystalline equilibrium phases of particle systems as well as their 
(zero-temperature) ground-state structures in low dimensions in which the interactions are 
characterized by steep repulsions and short-ranged attractions. 



The preponderance o 
in various dimensions 



of previous 

i3i,fly,3. 



rejvious investigations have focused on dense packings of spheres 
For congruent particles in three dimensions, the sphere 
is the only non-tiling particle for which the densest packing arrangements can be proved 
[lo| . It is only very recently that attention has turned to finding the maximal-density 
pac king arrangements of nonspherical particles in M^, including ellipsoids flll.[l^. tetrahedra 



1131,11 



17(1 . Very little is known about the densest packings of 



and superballs |16 
polyhedral particles jisl . 

The Platonic solids (mentioned in Plato's Timaeus) are convex polyhedra with faces 
composed of congruent convex regular polygons. There are exactly five such solids: the 
tetrahedron (PI), icosahedron (P2), dodecahedron (P3), octahedron (P4), and cube (P5) 
(see Fig. [1]) 19|. One major concern in this paper is the determination of the densest 
packings of each of the Platonic solids in three-dimensional Euclidean space M'^, except for 
the cube, which is the only Platonic solid that tiles space. 
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FIG. 1: (color online). The five Platonic solids: tetrahedron (PI), icosahedron (P2), dodecahedron 
(P3), octahedron (P4), and cube (P5). 

It is useful to highlight some basic geometrical properties of the Platonic solids that we 
will employ in subsequent sections of the paper. The dihedral angle 9 is the interior angle 
between any two face planes and is given by 

. 6 cos(7r/g) 

sm - = — — , (1) 

2 sin(7r/p) ^ ' 

where p is the number of sides of each face and q is the number of faces meeting at each 



vertex. Thus, ^ is 2 sin ^(1/V3), 2sin"^($/V3), 2 sin-i(<l>/V$2 + 1), 2sin^'\^/2j?,), and 
7r/2, for the tetrahedron, icosahedron, dodecahedron, octahedron, and cube, respectively, 
where $ = (1 + \/h)/2 is the golden ratio. Thus, since the dihedral angle for the cube is the 
only one that is a submultiple of 27r, the cube is the only Platonic solid that tiles space. We 
note in passing that in addition to the regular tessellation of space by cubes in the simple 
cubic lattice arrangement, there are an infinite number of other irregular tessellations of 



space by cubes [20| . Figure [2] shows a portion of a realization of a two-dimensional analog 
of such an irregular tessellation. 

Every polyhedron has a dual polyhedron with faces and vertices interchanged. The dual 
of each Platonic solid is another Platonic solid, and therefore they can be arranged into dual 
pairs: the tetrahedron is self-dual (i.e., its dual is another tetrahedron), the icosahedron and 
dodecahedron form a dual pair, and the octahedron and cube form a dual pair. 

An Archimedean solid is a highly symmetric, semi-regular convex polyhedron composed 
of two or more types of regular polygons meeting in identical vertices. There are thirteen 
Archimedean solids: truncated tetrahedron (Al), truncated icosahedron (A2), snub cube 
(A3), snub dodecahedron (A4), rhombicosidodecahdron (A5), truncated icosidodecahdron 
(A6), truncated cuboctahedron (A7), icosidodecahedron (A8), rhombicuboctahedron (A9), 



FIG. 2: (color online). A portion of a realization of an irregular tiling of the plane by squares. 




A1 A2 A3 A4 A5 




A6 A7 A8 A9 A10 




A11 A12 A13 

FIG. 3: The 13 Archimedean solids: truncated tetrahedron (Al), truncated icosahedron (A2), 
snub cube (A3), snub dodecahedron (A4), rhombicosidodecahdron (A5), truncated icosidodecah- 
dron (A6), truncated cuboctahedron (A7), icosidodecahedron (A8), rhombicuboctahedron (A9), 
truncated dodecahedron (AlO), cuboctahedron (All), truncated cube (A12), and truncated octa- 
hedron (A13). 
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truncated dodecahedron (AlO), cuboctahedron (All), truncated cube (A12) and truncated 
octahedron (A13) (see Fig. [3]). Note that the truncated octahedron is the only Archimedean 
sohd that tiles space. 

Another important observation is that the tetrahedron (PI) and the truncated tetrahe- 
dron (Al) are the only Platonic and Archimedean solids, respectively, that are not centrally 
symmetric. A particle is centrally symmetric if it has a center C that bisects every chord 
through C connecting any two boundary points of the particle, i.e., the center is a point of 
inversion symmetry. We will see that the central symmetry of the majority of the Platonic 
and Archimedean solids (P2 - P5, A2 - A13) distinguish their dense packing arrangements 
from those of the non-centrally symmetric ones (PI and Al) in a fundamental way. 

Some basic definitions concerning packings are given here. A saturated packing is one 
in which there is no space available to add another particle to the packing. A lattice A in 

is a subgroup consisting of the integer linear combinations of vectors that constitute a 



basis for M'^ [2l|. A lattice packing Pl is one in which the centroids of the nonoverlapping 
particles are located at the points of A, each oriented in the same direction. The set of 
lattice packings is a subset of all possible packings in M'^. In a lattice packing, the space M'^ 
can be geometrically divided into identical regions F called fundamental cells, each of which 
contains just the centroid of one particle. Thus, the density of a lattice packing is given by 



(2) 



^ Vol(F) ' 

where Vp is the volume of a ci-dimensional particle and Vol(F) is the volume of a fundamental 
cell. 

A more general notion than a lattice packing is a periodic packing. A periodic packing of 
congruent particles is obtained by placing a fixed nonoverlapping configuration of particles 
(where > 1) with arbitrary orientations in each fundamental cell of a lattice A. Thus, the 
packing is still periodic under translations by A, but the particles can occur anywhere in 
the chosen fundamental cell subject to the nonoverlap condition. The packing density of a 
periodic packing is given by 

- pvp, (3) 



Vol(F) 

where p = N/Vo\{F) is the number density, i.e., the number of particles per unit volume. 

The determination of the maximal-density arrangements of non-tiling polyhedral parti- 
cles is a notoriously difficult problem, especially since such extremal structures will generally 



be non-Bravais-lattice packings. Computer simulations that seek the maximal- density pack- 
ings can be an indispensable tool, especially if they can incorporate collective motions of the 
particles in order to obtain, in principle, the highest possible densities. However, the chal- 
lenge presented by polyhedral particles in is the non-smooth (i.e., nonanalytic) nature 
of the particle shape. In the case of smoothly-shaped particles, such as spheres, ellipsoids 
and superballs, one can construct analytic "overlap potential functions" for the particles 24 1 
and hence one can employ efficient collision- driven molecular dynamics (MP) growth pack- 



ing algorithms that inherently involve collective particle motions [1^, [161, [IT], |25|, |26| . The 
fact that analytic overlap potential functions cannot be constructed for polyhedral particles 
prevents us from using event-driven MD growth methods to study such systems. 

In this paper, we devise a novel optimization scheme, called the Adaptive Shrinking Cell 
(ASC), that can be applied to generate dense packings of polyhedra in M^. We employ it 
specifically to obtain the densest known packings of tetrahedra, icosahedra, dodecahedra, 
and octahedra with densities 0.823 . . . , 0.836 . . ., 0.904 . . ., and 0.947 . . ., respectively. The 



result for tetrahedra improves upon the density reported in our recent investigation 
Unlike the densest tetrahedral packing, which must not be a Bravais lattice packing 



27 



28|. 



ISj, the 



densest packings of the other non-tiling Platonic solids that we obtain are their corresponding 



densest lattice packings |22l. |23|. 



We also derive a simple upper bound on the maximal density of packings of congruent 
nonspherical particles, and apply it to Platonic solids, Archimedean solids, superballs and 
ellipsoids. We introduce the "asphericity" parameter 7 (ratio of the circumradius to inradius) 
to show that when 7 is sufficiently small, the upper bounds are relatively tight and thus close 
to the corresponding densities of the optimal lattice packings of octahedra, dodecahedra and 
icosahedra as well as of the majority of the Archimedean solids with central symmetry. 

Our simulation results as well as theoretical arguments lead us to conjecture that the 
densest packings of Platonic and Archimedean solids with central symmetry are given by 
their corresponding densest lattice packings. This can be regarded to be the analog of 
Kepler's sphere conjecture for these solids. The truncated tetrahedron is the only non- 
centrally symmetric Archimedean solid, the densest known packing of which is a non-lattice 
packing with density as high as 23/24 = 0.958333.... Our work also suggests that the 
optimal packings of superballs are their corresponding densest lattice packings. 

In a recent letter 27j, we briefly reported the densest known packings of the non-tiling 
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Platonic solids obtained using the ASC algorithm and proposed the aforementioned conjec- 
ture concerning the optimal packings of the Platonic and Archimedean solids. In this paper, 
we expand on theoretical and computational details and report additional new results. In 
particular, we provide comprehensive details about the ASC scheme and the simulation re- 
sults (Sec. II), including a discussion about the various initial configurations we used for the 
ASC algorithm (Sec. III]. Moreover, we have improved on the highest tetrahedral packing 



density reported in Ref. 



27( 1 (i.e., from 0.782 ... to 0.823 . . .), by exploring a broad range of 



dynamical parameters for the algorithm and initial configurations. Certain pair statistics of 
the densest known packing of tetrahedra (e.g., the contact number, the centroidal correla- 
tion function and the face- normal correlation function) are given (Sec. III). It is noteworthy 
that the densest tetrahedral packing is a non-Bravais structure with a complex periodic cell 
and possesses no long-range order. 

The initial configurations for the icosahedral, dodecahedral and octahedral packings are 
described and the numerical challenges in producing dense packings of such polyhedral 
particles are discussed (Sec. III). In addition, a detailed derivation of the upper bound and 
tables containing the geometrical characteristics of the Platonic and Archimedean solids as 
well as their upper bound values are given (Sec. IV). The upper bound is also applied 
to superballs and ellipsoids (which was not done in Ref. l^]). Moreover, we provide the 
major elements of a possible proof of our conjecture (Sec. V). We also discuss how our 
conjecture could be generalized to other centrally symmetric polyhedral particles such as 
prisms and anti-prisms as well as high- dimensional analogs of the Platonic solids. Our 
work also naturally leads to another conjecture reported here for the first time, namely, the 
optimal packing of any convex, congruent polyhedron without central symmetry is generally 
not a (Bravais) lattice packing (Sec. V). 

Furthermore, we discuss the possible applications and generalizations of the ASC scheme 
to predict the crystal structures of polyhedral nanoparticles and to the study of random 
packings of hard polyhedra. Finally, we collect in appendices basic packing characteristics 
of various optimal lattice and non-lattice packings of polyhedra (including lattice vectors) 
that have been scattered throughout the literature and provide lattice vectors and other 
characteristics of the densest known packings of tetrahedra (obtained here) and truncated 
tetrahedra. 
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II. ADAPTIVE SHRINKING CELL (ASC) OPTIMIZATION SCHEME 



We formulate the problem of generating dense packings of nonoverlapping, non-tiling 
polyhedra within a fundamental cell in subject to periodic boundary conditions as an 
optimization problem. In particular, the objective function is taken to be the negative of the 
packing density 0. Starting from an initial unsaturated packing configuration of particles of 
fixed size in the fundamental cell, the positions and orientations of the polyhedra are design 
variables for the optimization. Importantly, we also allow the boundary of the fundamental 
cell to deform macroscopically as well as compress or expand (while keeping the particles 
fixed in size) such that there is a net compression (increase of the density of the packing) in 



the final state 



291 ]. Thus, the deformation and compression/expansion of the cell boundary. 



which we call the adaptive fundamental cell, are also design variables. We are not aware of 
any packing algorithm that employs both a sequential search of the configurational space 
of the particles and the space of lattices via an adaptive fundamental cell that shrinks 
on average to obtain dense packings. We will call this optimization scheme the Adaptive 
Shrinking Cell (ASC). 

We will see that the ASC optimization scheme allows for some desired collective motions 
of the particles to find the optimal lattice for the periodic cell. This is to be contrasted with 
previous treatments that use a fixed shape for the fundamental cell, which may or may not 
be the optimal shape. Figure H] illustrates a simple sequence of configuration changes for a 
four-particle system (which is explained in more detail in Sec. II. A). By efficiently exploring 
the design- variable space (DVS), which consists of the particle configurational space as well 
as the space of lattices via an adaptive fundamental cell, the ASC scheme enables one to 
find a point in the DVS in the neighborhood of the starting point that has a higher packing 
density than the initial density. The process is continued until the deepest minimum of the 
objective function (a maximum of packing density) is obtained, which could be either a local 
or global optimum, depending on the particle shapes. 

The ASC optimization problem could be solved using various techniques, depending on 
the shapes of the particles. For example, for spheres, one can linearize the objective function 
and the constraints, and we have found that linear programming techniques can efficiently 
produce optimal solutions for such ASC problems 30|]. However, for hard polyhedra the 



non-overlapping conditions given by the separation axis theorem (discussed in detail below) 
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(a) 



(b) 



(c) 



FIG. 4: (color online). Sequential changes of the packing configuration due to the design vari- 
ables in the ASC algorithm, (a) An initial configuration of four particles, (b) A trial move of a 
randomly selected particle that is rejected because it overlaps another particle, (c) A trial move 
that is accepted, which results in a deformation and compression (small in magnitude) changing 
the fundamental cell shape and size as well as the relative distances between the particles. The 
large fundamental cell is divided into smaller sub-cells in order to implement the "cell method" 
discussed in Sec. II. D. 

involve at least quartic inequalities, which makes it inefficient to solve even using nonlinear 
programming methods. 

Here, for polyhedral particles, we solve the ASC optimization problem using a stochastic 
procedure, i.e., Monte Carlo (MC) method with a Metropolis acceptance rule for trial moves 
to search the DVS efficiently. However, it is important to distinguish our procedure that 
incorporates deformation and compression/expansion of the fundamental cell (i.e., t he sp ace 
of lattices) as design variables from previous MC hard-particle packing algorithms 3l|. In 
standard MC simulations, arbitrarily selected individual particle is given random displace- 
ment or rotation. This sequential movement method is not able to account for any collective 
motions of the particles, which is crucial to increasing the packing density, as pointed out in 
Sec. I. In our procedure, the deformation/compression/expansion of the boundary at least 
in part allows for collective particle motions in a direction leading to higher packing den- 
sity. Moreover, it is the overall compression of the fundamental cell that causes the packing 
density to increase, not the growth of the particles as in most MD and MC hard-particle 
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packing algorithms [25|, l26|, |3l| . 



At first glance, one might surmise that an algorithm that employs particle growth with 
an adaptive non-shrinking fundamental cell is equivalent to our choice of fixing the particle 
size while allowing the cell to shrink on average. The former is computationally less efficient 
than the latter for polyhedral particles. Specifically, growth of polyhedral particles requires 
manipulating the coordinates of the vertices of each particle, and thus involves at least dNn^ 
numerical operations, where d is the spatial dimension, is the total number of particles in 
the system and n„ is the number of vertices per particle. It is much more computationally 
expensive to use growing/shrinking particles as trial moves for the optimization scheme, 
especially when the number of particles is large, compared to our adaptive fundamental cell 
approach, which only requires manipulating d{d+ l)/2 strain components. 

In the ASC scheme, the macroscopic deformation and compression/expansion of the fun- 
damental cell of the lattice is completely specified by a strain tensor. Since we only consider 
small deformations, linear strain analysis can be applied here. Starting from an initial con- 
figuration, a trial configuration can be generated by moving (translating and rotating) a ran- 
domly chosen particle or by a random macroscopic deformation and compression/expansion 
of the fundamental cell. If any two particles overlap, the trial configuration is rejected; 
otherwise, if the fundamental cell shrinks in size (which makes the density higher), the 
trial configuration is accepted. On the other hand, if the cell expands in size, the trial 
configuration is accepted with a specified probability pacc, which decreases as increases 
and approaches zero when the jamming limit 321] (i.e., locally maximal dense packing) is 
reached. In particular, we find Pacc-, with an initial value pacc ~ 0.35, decreasing as a power 
law with exponent equal to -1 works well for most systems that we studied. The particle 
motion is equally likely to be a translation or a rotation. The ratio of the number of particle 
motions to the number of boundary trial moves should be greater than unity (especially 
towards the end of the simulation), since compressing a dense packing could result in many 
overlaps between the particles. Depending on the initial density, the magnitudes of the 
particle motions and strain components (e.g., the dynamical parameters) need to be chosen 
carefully to avoid the system getting stuck in some shallow local minimum. The parame- 
ters should also be adjusted accordingly as the packing density increases, especially towards 
the jamming point. The number of total Monte-Carlo moves per particle is of the order of 
5 X 10^ 
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In the ensuing subsections, we describe in detail how we implement particle motions as 
well as the deformation and compression/expansion of the fundamental cell, and precise 
check for interparticle overlaps using the separation axis theorem 331]. We also discuss the 
cell method and nearest-neighbor lists that are employed to speed up the simulations. 



A. Particle Motions 

The fact that a polyhedron is the convex hull of its vertices makes the set of vertices a 
useful geometrical representation of a such particle. It is convenient to choose the origin 
of the local coordinate system for the vertices to be the centroid of the polyhedron. Other 
important geometrical properties of the polyhedron, such as its faces and edges, can be 
represented as certain subsets of the vertices. 

Let the Euclidean position of the centroid of the jth particle be x|' {j = 1,..,N). A 
translational motion of the particle centroid can be obtained by generating a randomly 
oriented displacement Ax;^ with small magnitude (10~^ — 10"^ of the characteristic length 



of the particle), i.e., 



xf = x; + Ax; 



(4) 



A rotational motion can be generated by rotating the particle (all of its vertices) along a 
randomly chosen axis (passing through its centroid) by a random small angle 6. Let the 
vector Vj originating at the centroid denote the vertex i. We have 



R ■ v,^ + 



(5) 



where and v| are the components of Vj perpendicular and parallel to the rotation axis, 
respectively, and R is the rotation matrix, i.e.. 



R 



cos 6 — sin 6 
sin 6 cos 6 
1 



(6) 



If this motion does not result in the overlap with another particle, the trial move is accepted; 
otherwise it is rejected. 
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B. Adaptive Fundamental Cell 



For a lattice-based periodic packing, the fundamental cell is specified by the lattice vectors 
aj {i = 1,2,3). Recall that the Euclidean coordinates of the particle centroids are x|' 
(j = 1, .., A^). The relative coordinates of the centroids with respect to the lattice vectors 
are given by 



r . X- 



^ (7) 

where T = [ai, a2, as]. 

The adaptive fundamental cell allows for a small strain of the fundamental cell, including 
both volume and shape changes, which is represented by a symmetric strain tensor e, i.e.. 



AT = £ ■ r, 

where 

Cll ^12 ^13 
^21 £22 £23 

£31 £32 £33 

and the new fundamental cell (lattice vectors) are given by 

f = r + AT. 

Substituting the above equation into Eq. (jl]), we have 



(8) 



(9) 



r-x 



x^ + AT ■ x^ 



(10) 



(11) 



Thus, the strain of the fundamental cell corresponds to non-trivial collective motions of 
the particle centroids (see Fig. H]). In general, the translational motions of the particles con- 
tain the contributions from a random independent part (given by Eq. (jlj)) and the collective 
motion imposed by the adaptive fundamental cell. It is this collective motion that enables 
the algorithm to explore the configuration space more efficiently and to produce highly dense 
packings. 



C. Checking Overlaps and the Separation Axis Theorem 

Hard polyhedron particles, unlike spheres, ellipsoids and superballs, do not possess simple 
overlap potential functions. Thus, the check of overlapping for such particles requires other 
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techniques. Two convex objects are separated in space if and only if there exist an axis, on 
which the hne segments defined by projections of the two objects do not overlap (see Fig. [5]). 
This statement is usually referred to as the separation axis theorem (SAT) 




FIG. 5: (color online). Two nonoverlapping particles and one of their separation axes. 

For convex polyhedra, the above theorem has a simpler version. Since a polyhedron 
is completely defined by its vertices, the line segment on the axis is defined by the most 
separated two projections of the vertices. Moreover, the axis is either perpendicular to one 
of the faces or perpendicular to a pair of edges, one from each polyhedron. This reduces the 
number of axes that need to be checked from infinity to [E{E — l)/2 + 2F], where E and F 
is the number of edges and faces of the polyhedron, respectively. Here we employ the SAT 
to check for interparticle overlaps numerically up to the highest machine precision. 

For polyhedra whose circumscribed and inscribed spheres are well defined, two particles 
are guaranteed to overlap if the centroidal separation is smaller than the inscribed diameter 
and guaranteed not to overlap if the centroidal separation is larger than the circumscribed 
diameter. This pre-check dramatically speeds up the simulations starting from configurations 
with low densities. 



D. Cell Method and Near-Neighbor List 



The well-developed cell method 25|, |26| for particle-system simulation is used here to 
speed up the process. However, the simulation box (fundamental cell) will not remain a 
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FIG. 6: (color online). The particles (black) in the central fundamental cell (with red bound- 
ary) and their images (gray). The distances between the particles in the central cell and all the 
corresponding images need to be checked. 



cube during the simulation and several conventional techniques developed for cubic box 
need modifications. In particular, to obtain the minimum image distances, all surrounding 
boxes need to be checked explicitly (see Fig. [6]). With the help of the relative-coordinate 
representation of the particle centroids, the separation vector of a pair of particles is given 



where L is a vector with the values of components to be 1, or -1. Note L can also be 
considered as the index of the box that the image particle is in. 

The cells are taken to be the same shapes as the simulation box. Partitions of the particle 
centroids into the cells are also convenient for relative coordinate representations, i.e., the 
index of the cell (represented by a vector) C for particle i is given by 



where [X] gives the smallest integer part of X and Nc is the number of cells. Note that 
the boundary deformations and rotations of particles do not affect the cells in which the 
particles are situated; only translations may cause transitions between the different cells. 



by 



ds = T- (xf - + L) 



(12) 



C=[iVcxf], 



(13) 
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For dense packings, we ma^ 



efficiency of the algorithm 



ce use of a near-neighbor hst (NNL) to further improve the 



25 



261 . In particular, when the packing density is high, each 



particle is "trapped" in a "cage" formed by its near neighbors. These near-neighbor config- 
urations are practically stationary, i.e., the particles undergo very small "jiggling" motions. 
Thus, to check for overlapping between the particles, one only needs to consider a particles 
nearest neighbors. If a large magnitude translation is made, the NNL needs to be set up 
again. 



III. APPLICATION OF THE ASC SCHEME TO THE PLATONIC SOLIDS AND 
RESULTS 

Here we apply the ASC scheme to obtain the densest known packings of the non-tiling 
Platonic solids. Due to its lack of central symmetry, the tetrahedron presents the greatest 
challenge for the numerical solution procedure of the ASC scheme because of its tendency 
to get stuck in local (density) minima, which is a consequence of the associated "rugged" 
energy or, more precisely, "density" landscape (e.g., the packing density as function of the 
centroidal positions and orientations of all of the particles in the fundamental cell). Hence, 
the choice of initial configurations becomes crucial in getting dense tetrahedral packings. By 
contrast, the central symmetry of the octahedron, dodecahedron and icosahedron results in 
density landscapes for the numerical procedure that are appreciably less rugged than that 
for tetrahedral packings. 



A. Tetrahedra 



The determination of the densest packings of regular tetrahedra is part of the 18th prob- 
lem of Hubert's famous set of problems. It is of interest to note that the densest (Bravais) 
lattice packing of tetrahedra (which requires all of the tetrahedra to have the same orienta- 
tions) has the relatively low density 0max = 18/49 = 0.367 . . . and each tetrahedron touches 
14 others 3^. Recently, Conway and Torquato showed that the densest packings of regular 
tetrahedra must not be a Bravais lattice packings, and found packings with density as large 
as ~ 0.72 One such packing is based upon the filling of "imaginary" icosahedra with 
the densest arrangement of 20 tetrahedra and then employing the densest lattice packing 
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"Welsh" packings with density (j) 

et al. 1^ experimentally generated jammed disordered packings of such dice with 



of icosahedra. A slightly higher density was achieved by a perturbation of the so-called 

17/24 ^ 0.708 Q. Using "tetrahedral" dice, Chaikin 

a 0.75. 

Even though the dice are not perfect tetrahedra (because the vertices and edges are slightly 
rounded), these experimental results suggested that the densest packings of tetrahedra could 
exceed the highest densities reported by Conway and Torquato. Indeed, Chen [l5| has re- 
cently discovered a periodic packing of tetrahedra with (p = 0.7786. . . js^. We will call 
this the "wagon-wheels" packing because the basic subunits consist of two orthogonally in- 



tersecting "wagon" wheels. A "wagon whee 
around a common edge (see Fig. la of Ref. 



consists of five contacting tetrahedra packed 



TABLE I: Some packing characteristics of certain known tetrahedral packings. Here (j) is the 
packing density, is the number of particles in the fundamental cell, and Z is the average contact 
number per particle. 



Name 


<P 


Locally Jammed 


N 


Z 


Optimal Lattice [34] 


18/49 « 0.367346 


Yes 


1 


14 


Uniform [131 


2/3 w 0.666666 


Yes 


2 


10 


Welsh [13] 


17/24 0.708333 


No 


34 


25.9 


Icosahedral [131 


0.716559 


No 


20 


20.6 


Wagon Wheels [15] 


0.778615 


Yes 


18 


7.1 



1. Initial Conditions 

Because we will use all of the aforementioned tetrahedral packing arrangements as ini- 
tial conditions in our numerical solution procedure for the ASC optimization problem, we 
describe them now in a bit more detail. Table [T] summarizes some of their packing charac- 
teristics, including the packing density 0, the number of particles in the fundamental cell, 
the average number of contacting particles per particle Z, and whether the packing is locally 
jammed. A packing is locally jammed if each particle is locally trapped by its neighbors, 
i.e., it cannot be translated or rotated while fixing the positions and orientations of all the 



other particles 32|, |36| 
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In the optimal lattice packing, each tetrahedron touches 14 others via edges and vertices. 
In the uniform packing, the tetrahedra are "locked" in a planer layer, and each particle 
has eight in-plane partial face-to-face contacts and two edge-to-edge contacts (contributed 
by the two layers above and below), and therefore each particle contacts 10 others. In the 
icosahedral packing, the tetrahedra have a greater degree of freedom to move, and indeed it 
is not locally jammed. The 20 tetrahedra inside the imaginary icosahedron are required to 
meet at its centroid, and so each tetrahedron has 20 vertex-to-vertex contacts. Placing the 
imaginary icosahedra on the sites of its optimal lattice results in a tetrahedral packing with 12 
partial face-to-face contacts. Thus, there are 8 tetrahedra that have 20 contacts per particle 
and 12 tetrahedra with 21 contacts per particle, and therefore Z = (8 x 20 -|- 12 x 21)/20 = 
20.6. 

The Welsh packing is based on the "primitive Welsh" tessellation of M'^ into truncated 
large tetrahedra and small regular tetrahedra 1^ • Interestingly, the "primitive Welsh" con- 
figuration is closely related to the uniform packing, i.e., it can be constructed by replicating 
a periodic packing with two truncated large tetrahedra and two small regular tetrahedra 
in the fundamental cell. Each truncated large tetrahedron can be divided into 12 "Welsh 
lows", and 4 "Welsh medials," following the notation in Ref. 1^. The Welsh medials and 
Welsh lows are tetrahedra with lower symmetry than regular tetrahedra, which are referred 
to as the "Welsh highs." To construct the Welsh packing, we insert Welsh highs (regular 
tetrahedra) into the Welsh low and medial regions. We will refer to the regular tetrahedra 
that are placed in the Welsh low, medial and high regions as lows^ medials and highs, re- 
spectively. In the Welsh packing, each high makes 4 face-to-face contacts with 4 medials. 
Each medial contacts 24 lows and 3 other medials at the centroid of the large truncated 
tetrahedron; moreover, each medial makes a face-to-face contact with one high. Thus, each 
medial has 28 contacts. Each low contacts 23 other lows and 4 medials at the centroid of 
the large truncated tetrahedron. Thus, in total each low has 27 contacts. Each truncated 
large tetrahedron has 12 lows and 4 medials. Thus, the number of basis particles in the 
fundamental cell = 34 and Z = (12 x 27 + 4 x 28 + 1 x 4)/17 = 440/17 = 25.882352 . . .. 

In the fundamental cell of the wagon- wheels packing, there are = 18 particles, forming 
two clusters, each of which includes two "wheels" entangled orthogonally to each other. 
Each cluster has a central connection particle, an "upper" wheel and a "lower" wheel. 
The central particle contacts the other 8 tetrahedra through two of its edges. Each of 
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the other 8 tetrahedra has 4 edge contacts with particles in its own cluster and 3 partial 
face-to-face contacts with particles from other clusters. Thus the average contact number 
Z = (1 X 8 + 8 X 7)/9 = 64/9 = 7.111111 . . .. The coordinates of all of the 18 tetrahedra in 
the fundamental cell can be found on the authors' website 37 1. 



2. Dense Packings 

We employ our algorithm to obtain dense packings of tetrahedra using initial configu- 
rations that are based upon the known packings given in Table [T] as well as certain dilute 
packings with carefully chosen fundamental cells and the number of particles. Note that 
these dilute packings were not used as initial configurations to obtain the results reported 

n 

in Ref. [27|. A range of initial densities is used to yield the largest possible densities in the 
final states, as summarized in Table [TTl 

For initial conditions using the known packings in Table [H the fundamental cell is isotropi- 
cally expanded with the relative coordinates of the tetrahedra fixed so that the initial packing 
is locally unjammed with a lower density (pint- The highest density of 0.782 obtained from 
this subset of initial conditions started from the "wagon-wheels" packing, which was the 



value reported in Ref. [27 1. 



For the aforementioned dilute packings, a variety of fundamental cell shapes spanning 
from that for the simple cubic lattice to that of the hexagonal close packing as well as a 
wide range of particle numbers spanning from 70 to 350 are explored to generate the dense 
packings. Specifically, the densest known tetrahedral packing is obtained from a dilute initial 
packing with a rhombical fundamental cell similar to that of the hexagonal close packing 
and 314 particles (see Table HB . The particles are originally placed in the fundamental cell 
randomly and then the system is sufficiently equilibrated with the fundamental cell fixed. 
The final packing configuration has a density of about 0.823 and is shown in Fig. [7] from 
two different viewpoints. The optimized lattice vectors for the densest packing and other 
characteristics are given in Appendix B. We see that the packing lacks long-range order 
and is composed of '"clusters" of distorted wagon wheels and individual tetrahedra. Thus, 
we call this dense arrangement the "disordered wagon wheels" packing. We will comment 
further about the significance of achieving this remarkably high density with a disordered 
packing in the Discussion and Conclusions (Sec. V). 
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TABLE II: Dense tetrahedral packings generated from our algorithm using initial configurations 
based upon the packings given Tableland the dilute packing as described in the text. Here (pint 
is the initial packing density, (j) is the final packing density, is the number of tetrahedra in the 



fundamental cell 



381], and Z is the average contact number per particle. 



Initial Packing 



(pint 



Locally Jammed N Z 



Optimal Lattice 
Uniform 
Welsh 
Icosahedral 
Wagon Wheels 



0.25 - 0.3 0.695407 

0.45 - 0.55 0.665384 

0.45 - 0.6 0.752092 

0.45 - 0.6 0.744532 

0.55 - 0.65 0.782021 



Disordered Wagon Wheels 0.005 - 0.01 0.822637 



Yes 
Yes 
Yes 
Yes 
Yes 
Yes 



27 6.6 

54 9.8 

34 7.4 

20 7.1 

72 7.6 

314 7.4 




(a) (b) 

FIG. 7: The densest known packing of tetrahedra with 314 particles in the fundamental cell. We 
call this the "disordered wagon-wheels" packing, (a) Viewed from the side, (b) Viewed from the 
top. It is apparent that the packing lacks long-range order and consists of "clusters" of distorted 
wagon wheels. 

3. Statistics of the Densest Tetrahedral Packing 

Contact Numbers: The contact numbers of the densest tetrahedral packing as well as 
those obtained using the other five initial configurations (Table [T]) are given in Table [Tll To 
determine contacting neighbors, it is crucial to find the interparticle gaps first. In particular, 
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each of the equilateral triangular faces of the tetrahedra is discretized into M equal-sized 
smaller triangles and the distances between the vertices of the small triangles of the neigh- 
boring tetrahedra are computed, the minimum of which is used as the gap value between 
the two corresponding tetrahedra. For our packings, when M > 200, the values of interpar- 
ticle gaps obtained by this approximation scheme become stable and converge to the true 
gap values. Like any packing generated via computer simulations, the particles never form 
perfect contacts and hence contacting neighbors must be determined by setting a tolerance 
T for the interparticle gaps. Here we choose T to be equal to the mean value of the gaps 
(of order 10^^ to 10~^ of the edge length) such that any gap less than T is associated with 
a contact. This procedure yields the contact numbers given in Table [TTl In particular, we 
find that in the disordered wagon wheels packing, there are approximately 6.3 face-to-face 
or partial face-to- face contacts and 1.1 edge-to-edge contacts. 

Face-Normal Correlation Function: An important statistical descriptor for packings of 
non-spherical particles is the particle orientation correlation function, which measures the 
extent to which a particle's orientation affects the orientation of another particle at a different 
position. For highly symmetric particle shapes, it is reasonable to focus on those pair 
configurations in which the particles are in the same orientation, since particle alignment 
is associated with dense configurations and phase transitions that may occur. However, 
because tetrahedra lack central symmetry, we found that face-to-face contacts are favored 
by the dense packings, instead of face-to-vertex contacts, which are necessarily associated 
with aligning tetrahedral configurations, such as in the case of the optimal lattice packing. 




12 3 4 



FIG. 8: The face-normal correlation function CpNir) of the densest known packing of tetrahedra 
obtained here, where rj„ is the radius of the insphere of a tetrahedron. 
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Thus, because our interest is in dense tetrahedral packings, we determine the "face- 
normal" correlation function CpNir), which we define as the average of the largest negative 
value of the inner product of two face normals of a pair of tetrahedra separated by a distance 
r. This quantity is plotted in Fig. [HI The strong positive peak values, ranging from r = 
2rin = c/q/v^ (where rj„ is the radius of the insphere of the tetrahedron and do is the 
edge length of the particle) to r slightly greater than 2rm, indicate a large number of face- 
to-face contacts between neighboring particles, implying strong short-range orientational 
correlations. The fact that CpNif) plateaus to its maximum or near maximum value, for 
the distance interval r/{2rin) G (1, 2.2), indicates that there are partial face-to-face contacts 
in which particles slide relative to one another such that the distance between the centroids 
can increase while the face-normal inner product remains the same. The "valleys" of CpNif) 
with relatively small magnitudes are manifestations of weak particle alignments (alignment 
in the same direction), which are local configurations that are necessary "costs" to achieve 
the largest face-to-face contact numbers on average. 

Centroidal Radial Distribution Function: Another statistical descriptor of the structure 
is the centroidal radial distribution function g2ir). In particular, g2{r)r'^dr is proportional 
to the conditional probability that a particle centroid is found in a spherical shell with 
thickness dr at a radial distance r from another particle centroid at the origin. It is well 
established that when there is no long-range order in the system, g2 decays to unity very fast. 
Figure [9] shows the centroidal radial distribution of the densest known tetrahedral packing. 
We see that g2{j') decays to unity after a few oscillations, indicating that the packing lacks 
long-range order. 

B. Octahedra 

To obtain dense packings of octahedra, we use a wide range of initial configurations. This 
includes unsaturated packings in which the particles are randomly oriented and positioned 
with densities spanning the range from 0.2 to 0.35 as well as unsaturated simple cubic (SC) 
and face-centered cubic (FCC) lattice packings with densities from 0.1 to 0.125, and the 
optimal lattice packings with densities from 0.3 to 0.55 are used. We found that, using small 
compression rates for the random initial configurations and moderate compression rates for 
the various lattice configurations along with a sufficient number of particle displacements 
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FIG. 9: The centroidal radial distribution function 52 (^) of the densest known packing of tetrahedra, 
where rj„ is the radius of the insphere of a tetrahedron. 

and rotations, resulted in final configurations with densities larger than 0.93 that are very 
close in structure and density to the optimal lattice packing. More specifically, starting 
from the unsaturated SC lattice (with density 0.1) and optimal lattice (with density 0.55) 
packings, final packings with densities very slightly larger than the value 0.947003 . . . can be 
obtained, which are extremely close in structure and density to the optimal lattice packing 



{'Pmax — 0.947368. . .) [22I]. A previous study involving an event-driven molecular dynamics 
growth algorithm for the case of octahedral- like superballs led to the same the densest lattice 
packing, which lends further credence to the fact that this lattice is indeed optimal among all 
packings [l^] • Each octahedron in the optimal lattice packing, depicted in Fig. [T0|, contacts 
14 others and its lattice vectors are given in Appendix A. 




FIG. 10: A portion of the optimal lattice packing of octahedra. 
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C. Icosahedra 



As in the case of octahedra, we use a variety of initial conditions to get dense packings 
of icosahedra. This includes unsaturated packings in which the icosahedra are randomly 
oriented and positioned with densities spanning the range from 0.2 to 0.3 as well as various 
unsaturated lattice packings, such as the BCC, FCC and the densest lattice packing, with 
densities spanning the range from 0.3 to 0.65. We found that by starting from a low-density 
configuration, with a small enough compression rate and sufficient number of particle moves, 
one can always obtain a final configuration that is very close to the optimal lattice packing 
with a density ^ 0.83. An initial configuration of an unsaturated optimal lattice packing 
configuration with a density 0.65 gives a final packing with density 0.836315. . ., which is 
extremely close in structure and density to the optimal lattice packing {4>^ax = 0.836357 . . .) 



231]. Each icosahedron of the optimal lattice packing, depicted in Fig. [TTl contacts 12 others 



and its lattice vectors are given in Appendix A. 



FIG. 11: A portion of the optimal lattice packing of icosahedra. 



D. Dodecahedra 

Unsaturated random packings with densities ranging from 0.15 to 0.3 as well as unsatu- 
rated simple cubic and the optimal lattice packings with densities spanning the range from 
0.3 to 0.6 are employed as initial configurations to generate dense packings of dodecahedra. 
We found that it is algorithmically more difficult to avoid local density maxima for dodecahe- 
dra than octahedron and icosahedron packings discussed above. For example, starting from 
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random configurations, even with sufficiently small compression rate and a large number of 
particle moves, we can only achieve apparently jammed final packings with densities in the 
range from 0.83 to 0.85. When we use an unsaturated optimal lattice packing with density 
0.72 as an initial condition, we can generate a final packing with cf) = 0.904002 . . ., which is 
relatively close in structure and density to the optimal lattice packing {(t>f^ax — 0.904508 . . .) 
23|. Each dodecahedron of the optimal lattice packing, depicted in Fig. [121 contacts 12 



others and its lattice vectors are given in Appendix A. 




FIG. 12: A portion of the optimal lattice packing of dodecahedra. 

We note that the tendency for dodecahedral packings to get stuck in local-density maxima 
in our algorithm is due to the fact that in some sense the dodecahedron is a shape that is 
intermediate between the octahedron and icosahedron. Specifically, for octahedral packings, 
the octahedral symmetry of the particles facilitates the formation of nematic phases at 
relatively high densities, which can then be easily compressed into a dense crystalline phase 
with adaptive fundamental cells. An icosahedron is highly isotropic and possesses a large 
number of faces. Thus, at low densities, the icosahderal packing behaves like a hard-sphere 
fluid, and it is only near the jamming point that the asphericity of the icosahedron begins 
to play an important role. Again, by allowing the fundamental cell to adapt in the case of 
icosahedra, the optimal lattice that best accommodates the packing can be easily identified. 
However, for our algorithm, the dodecahedral packing behaves differently from either the 
octahedral or icosahedral packings. In particular, while dodecahedra are more isotropic 
than octahedra, which strongly suppresses the formation of nematic phases, they have fewer 
but larger faces than icosahedra, which favors the formation of face-to-face contacts that 
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ends up jamming the dodecahedral packing at lower densities than either the octahedral or 
icosahedral packings. 



IV. UPPER BOUND ON THE MAXIMAL DENSITY OF PACKINGS OF NON- 
SPHERICAL PARTICLES 

Here we derive a simple upper bound on the maximal density (pmax of a packing of 
congruent nonspherical particles of volume Vp in any Euclidean space dimension d. We will 
see that this bound will aid in our analysis of the optimality of not only the Platonic and 
Archimedean solids but also superballs. Let (p^ax be the maximal density of a d-dimensional 
packing of congruent spheres and let Vg represent the volume of the largest sphere than can 
be inscribed in the nonspherical particle. 

Lemma: The maximal density of a packing of congruent nonspherical particles is bounded 



from above according to the following bound 

< min 



■'max 



^ 6^ 1 

^maxi 



(14) 



where min[x, y\ denotes the minimum of x and y. 

The proof is straightforward. The maximal packing density (pmax can be expressed in 
terms of the maximal number density Pmax via the relation 



■'max 



pmax • 

(15) 



If we inscribe within each nonspherical particle of the packing the largest possible sphere, it 
is clear that 

pmaxVs < (Piax, (16) 

and therefore combination of the last two equations yields the upper bound of the Lemma. 
Remark: The upper bound ((Hj) will yield a reasonably tight bound for packings of non- 
spherical particles provided that the asphericity ^ of the particle is not large. Here we define 
the asphericity as the following ratio: 

7 = — (17) 

'"in 

where Tout and rj„ are the circumradius and inradius of the circumsphere and insphere of the 
nonspherical particle. The circumsphere is the smallest sphere containing the particle. The 
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insphere is the largest sphere than can be inscribed in the particle. For a sphere, clearly the 
asphericity 7 = 1. Since upper bound f|T^ cannot be sharp (i.e., exact) for a nonspherical 
particle, any packing construction for a nonspherical particle whose density is close to the 
upper bound fll4p is nearly optimal, if not optimal. 

In the three-dimensional case, the upper bound (HM becomes 



^max < 0^ax = ^in 



TT 



(18) 



V18 

We now apply this upper bound to packings of the Platonic and Archimedean solids and 
compare the bounds to the densities of the corresponding densest lattice packings. Moreover, 
we apply the upper bounds to superball and ellipsoid packings. 



A. Platonic Solids 



Table IIIII compares the density of the densest lattice packings of the Platonic solids to the 
corresponding upper bounds on the maximal density for such packings. The large asphericity 
and lack of central symmetry of the tetrahedron is consistent with the large gap between the 
upper bound density and densest lattice packing density, and the fact that there are non- 
lattice packings with density appreciably greater than (p^^x- other hand, the central 
symmetry of the octahedron, dodecahedron and icosahedron and their associated relatively 
small asphericities explain the corresponding small differences between (p^^x ^max 
is consistent with our simulation findings that indicate that their optimal arrangements are 
their respective densest lattice packings. 



B. Archimedean Solids 



We also compute the upper bound ( 1181) for each of the 13 Archimedean solids and compare 



them to the densities of the corresponding densest lattice packings [23|, l3J] . (In Appendix 
A we provide the lattice vectors for the optimal lattice packings of the Archimedean solids). 
Table IIVI summarizes the upper bounds on the maximal density for such packings. Not 
surprisingly, the truncated tetrahedron (the only Archimedean solid that is not centrally 
symmetric) has a large asphericity, implying that there are denser non-lattice packings, as 
we explicitly identify in Section |Vl The central symmetry of the majority of the Platonic 
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TABL 



23 



ni: Comparison of the densities of the densest lattice packings for the Platonic solids 
to the corresponding upper-bound densities as obtained from (|18|) . Here Vp is volume 
of the polyhedron with unit edge length, and rout are the radii of the insphere and circumsphere 
of the polyhedron with unit edge length, respectively, 7 = Tout /Tin is the asphericity, (pmax is the 
density of the optimal lattice packing and <i)^ax is the upper bound (fTH]) . The numerical values are 
reported up to the sixth decimal place. The naming code used here is the same one used in Fig. 1. 



Name 


Vp 


Tin 


Tout 


7 


Ymax 


r'max 



Tetrahedron (PI) if if ^ 3 0.367346 1 

Icosahedron (P2) MsWs) sV^WTE Vio+2x/5 1258410 0.836357 0.893417 



Dodecahedron (P3) V250+ii0v^ v^+v^ 1.258410 0.904508 0.981162 

Octahedron (P4) ^ ^ ^ 1.732050 0.947368 1 

Cube (P5) 1 i ^ 1.732050 1 1 

and Archimedean solids and their associated relatively small asphericities explain the cor- 
responding small differences between 4>^ax (p^ax is consistent with our simulation 
findings that strongly indicate that their optimal arrangements are their respective densest 
lattice packings. 



C. Super balls 

The upper bound ( ITSl) in the case of superballs p/7| can be expressed analytically for 
all values of the deformation parameter p. A three-dimensional superball is a centrally 
symmetric body in occupying the region 

|x2pP+ |x3pP < 1, (19) 

where Xi {i = 1,2,3) are Cartesian coordinates and p > is the deformation parameter, 
which indicates to what extent the particle shape has deformed from that of a sphere {p = 1). 
A superball can possess two types of shape anisotropy: cubic-like shapes for 1 < p < 00, with 
p = 00 corresponding to the perfect cube, and octahedral-like shapes for < p < 1, with 
p = 1/2 corresponding to the perfect regular octahedron. In Ref. Il7|, event-driven molecular 



dynamics growth algorithms as well as theoretical arguments led to the conjecture that 
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TABLE IV: Comparison of the densities of the densest lattice packings for the Archimedean sohds 
to the corresponding upper-bound densities as obtained from (|18p . Except for the cubeoctahedron 
34l |. the densities of the densest lattice packings were obtained by Betke and Henk 



23( 1 ■ Here Vp 



is volume of the polyhedron with unit edge length, rj„ and rout are the radii of the insphere and 
circumsphere of the polyhedron with unit edge length, respectively, 7 = rout/ fin is the asphericity, 
'Pmax is the density of the optimal lattice packing, (/>max the upper bound ([18]) and t = [1 + 
(19 — 3\/33)^ + (19 + 3\/33)3]/3 ~ 1.83929. . . is the Tribonacci constant. The numerical values 
are reported up to the sixth decimal place. The naming code used here is the same one used in 
Fig. 3. 



Name 



fout 



7 



-'max 



Truncated Tetrahedron (Al) 
Truncated Icosahedron (A2) 
Snub Cube (A3) 
Snub Dodecahedron (A4) 
Rhombicicosidodecahedron (A5) 
Truncated Icosidodecahedron (A6) 
Truncated Cuboctahedron (A7) 
Icosidodecahedron (A8) 
Rhombicuboctahedron (A9) 
Truncated Dodecahedron (AlO) 
Cuboctahedron (All) 
Truncated Cube (A12) 
Truncated Octahedron (A13) 



23^2 
12 

125+43^5 
4 

3v^f^+4v^t+T 
3^2^ 

37.616654 

95 + 50^5 

60+29^5 
3 

12+10\/2 
3 

45+17^5 
6 

22 + 14^2 

5(99+47^5) 
12 

5^2 
3 

21+14^2 



4 



\/21+9V5 -y/sS+lsVS 
2V2 4 



i-1 



4{2-t) 

1.980915 
3.523154 
2.016403 

2 



3-t 



8^2 



5+2^5 
5 

l+2v^ 
2 

V25+llv^ 
2v^ 
V2 



I+V2 
2 

2 



4(2-t) 

2.155837 

■v/31+12V5 
2 

•v/ll+4V5 
2 

\/5+2V2 
2 

1+V5 
2 

a/i3+672 
2 

a/74+30V5 
4 

1 

\/7+4V2 



10 
2 



1.914854 
1.092945 
1.175999 
1.088303 
1.079258 
1.107392 
1.158941 
1.175570 
1.210737 
1.192598 
1.414213 
1.473625 
1.290994 



0.680921 
0.784987 
0.787699 
0.788640 
0.804708 
0.827213 
0.849373 
0.864720 
0.875805 
0.897787 
0.918367 
0.973747 
1 



1 

0.838563 
0.934921 
0.855474 
0.835964 
0.897316 
0.875805 
0.938002 
1 

0.973871 
1 
1 
1 



the densest packings of superballs for all convex shapes (1/2 < p < 00) are certain lattice 
packings, depending on the value of the deformation parameter p. 

For convex superballs in the octahedral regime (0.5 < p < 1), the upper bound (fT8|) is 
explicitly given by 

^l..-^r^^*^B(LJJi±l\s(Lj^\ (20) 



108p^ \2p 2p J \2p p 
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similarly, for superballs in the cubic regime (p > 1), the upper bound (fT8|) is given by 



max 



\2p 2p J \2p p 

where B{x, y) = r{x)r{y)/r{x + y) and r(x) is the Euler Gamma function. Here and in the 
subsequent paragraphs concerning ellipsoids, we only state the nontrivial part of the upper 
bound ( fTSl) . For p near the sphere point {p = 1), the densest lattice packings l3| have 
densities that lie relatively close to the corresponding upper-bound values. For example 
for p = 0.99, 0.98 and 0.97, = 0.745327 . . . , 0.750274 . . ., and 0.755325 . . ., respectively 
which is to be compared to 4>max = 0.740835 . . . , 0.741318 . . ., and 0.741940 . . ., respectively. 
For p = 1.01, 1.02 and 1.03, = 0.747834 . . . , 0.755084 . . ., and 0.762233 . . ., respectively, 
which is to be compared to </)^„^. = 0.741720 . . . , 0.742966 . . ., and 0.744218 . . ., respectively. 
The fact that the constructed densest lattice packings of superballs have density 
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relatively close to the upper bound value 0^ax strengthens the arguments made in Ref. 
that suggested that the optimal packings are in fact given by these arrangements. If the 
optimal packings around p = 1 are indeed the lattice packings, it would be surprising that 
for other values of p a continuous deformation of the superballs would result in a transition 
from lattice packings to denser non-lattice packings. 



D. Ellipsoids 

In the case of ellipsoids in which the the ratio of the three semiaxes are given by 1 : a : /3 
{a,P > 1), the upper bound (fT8|) becomes 

01, = a(3^. (22) 

V io 

For prolate spheroids /3 = 1, and the bound becomes 



max 



a^. (23) 



For oblate spheroids with a = (3 > 1, the bound is given by 



a 



max \/l^ 



(24) 



The upper bounds for ellipsoids generally do not work as well as those for the centrally 
symmetric Platonic and Archimedean solids, and superballs. This is due to the fact that the 
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three principle directions (axes) are not equivalent for an ellipsoid, i.e., it generally possesses 
different semiaxes. Note that the three equivalent principle (orthogonal) axes (directions) 
of a centrally symmetric particle are those directions that are two-fold rotational symmetry 
axes such that the distances from the particle centroids to the particle surfaces are equal. 
In addition, the asphericity of an ellipsoid increases linearly as the largest aspect ratio a. 
increases without limit. By contrast, the three principle directions are equivalent for the 
centrally-symmetric Platonic and Archimedean solids as well as for superballs, for which the 
asphericity is always bounded and close to unity (see Tables HIT] and [IV|) . Thus, we see that 
an asphericity value close to unity is a necessary condition to have a centrally symmetric 
particle in which the three principle directions are equivalent. 



V. DISCUSSION AND CONCLUSIONS 



We have formulated the problem of generating dense packings of nonoverlapping poly- 
hedra within an adaptive fundamental cell subject to periodic boundary conditions as an 
optimization problem, which we call the Adaptive Shrinking Cell (ASC) scheme. The pro- 
cedure allows both a sequential search of the configurational space of the particles and the 
space of lattices via an adaptive fundamental cell that shrinks on average to obtain dense 
packings. We have applied the ASC to generate the densest known packings of the Platonic 
solids. 

For tetrahedra, we find a packing with density ~ 0.823, which is a periodic (non- 
Bravais lattice) packing with a complex basis. Unlike the other Platonic solids, finding 
dense packings of tetrahedra with our algorithm requires having good initial configurations. 
The densest packing was found using 314 particles in a rhombical fundamental cell that 
is similar to that of the hexagonal close packing. As we stressed in our earlier work 
and continue to confirm in this paper, it is possible that denser packings of tetrahedra will 
involve increasingly larger numbers of particles in the fundamental cell. In fact, the higher 
density found here is realized by a larger periodic packing (314 particles per cell) than 
the one reported in Ref. (72 particles per cell) with density 0.782. It is apparent the 
obtained densest known tetrahedral packing is disordered in the sense that it possesses no 
long-range order, at least on the scale of the simulation box. This packing can be considered 
to be a disordered "mixture" of distorted wagon wheels and individual tetrahedra and it 
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is distinct from the one reported in Ref. [27| in both the arrangement of the wagon wheels 
and the extent of the distortion of the wagon wheels. Although we cannot rule out the 
possibility of the existence of denser packings involving a more ordered arrangement of 
the wagon wheels and individual tetrahedra, it is reasonable to expect that wagon wheels 
will be key building blocks in denser packings and should at least be slightly distorted to fill 
the interparticle gaps in wagon-wheel clusters more efficiently, which necessarily introduces a 
certain degree of disorder to the packing. If denser packings of tetrahedra must involve larger 
number of particles without long-range order, than it raises the amazing prospect that the 
densest packings of tetrahedra might be truly disordered, due to the geometrical frustration 
associated with the lack of central symmetry of a tetrahedron and that tetrahedra cannot 
tile space. This would be the first example of a maximally dense packing of congruent convex 
three-dimensional particles without long-range order. However, we cannot offer definitive 
conclusions about this possibility at this stage. It is clear that in future work in the search 
for denser packings, increasingly larger number of tetrahedra must be considered, which can 
only be studied using greater computational resources. 

Our simulation results and rigorous upper bounds strongly suggest that the optimal 
lattice packings of the centrally symmetric Platonic solids (octahedra, dodecahedra and 
icosahedra) are indeed the densest packings of these particles, especially since these arise 
rom a variety of initial "dilute" multi-particle configurations within the fundamental cell 
4ol |. It is noteworthy that a different simulation procedure (an event-driven molecular 
dynamics growth algorithm with an adaptive fundamental cell) has recently been used to 
demonstrate that the densest packings of octahedral-like superballs (which contains the 
perfect octahedron) are likely to be the optimal lattice packings 17|. Moreover, the fact that 
the optimal lattice packing densities of certain centrally symmetric nonspherical particles, 
such as the Archimedean solids with central symmetry and convex superballs, are relatively 
close to their upper bounds as well as other theoretical arguments given below suggest that 
the densest packings of these particles may also be given by their optimal lattice packings. 

It is crucial to stress that the nonspherical particles in this family do not deviate appre- 
ciably from a sphere, i.e., their corresponding asphericity 7 is always bounded and relatively 
close to unity, and, moreover, their three principle axes (directions) are equivalent. This is 
in contrast to the ellipsoid, which, although centrally symmetric, generally possesses three 
principle axes (directions) that are inequivalent and its asphericity can increase without limit 
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as the largest aspect ratio grows. These characteristics suggest that the densest elhpsoid 
arrangements are non-lattice packings, which indeed has been verified 1^ . 




Our simulation results and the ensuing theoretical arguments lead to the following con- 
jecture: 

Conjecture 1: The densest packings of the centrally symmetric Platonic solids are given by 
their corresponding optimal lattice packings. 

We now sketch what could be the major elements of a proof of this conjecture. In the case 
of each of the Platonic solids, face-to-face contacts are favored over vertex-to-face contacts to 
achieve higher packing densities, since the former allows the particle centroids to get closer 
to one another. Such contacting neighbor configurations around each particle reduce the 
volume of the corresponding convex hull joining the centroids of the contacting particles, 
and the fraction of space covered by the particles within this convex hull should be increased; 
see Fig. [13] for a two-dimensional illustration. 

To achieve the densest packing, the fraction of space covered by the particles within the 
convex hull should ideally be maximized and so should the number of face-to-face contacts 
per particle. Of course, this is only a local criterion that may not be consistent with the 
densest global packing. However, it will be seen that the equivalence of the three princi- 
ple directions of a centrally symmetric Platonic solid is crucial for this local optimization 
criterion to be consistent with the globally densest packing. 

It is noteworthy that orienting each of the particles in the packings of centrally symmetric 
Platonic solids enable a larger number of face-to-face contacts and thus allows the maximal 
fraction of space covered by the particles within the convex hull. For example, a particle 
with F faces, possesses F/2 families of axes that go through the centroid of the particle and 
intersect the centrally-symmetric face pairs such that the particles (in the same orientation) 
with their centroids arranged on these axes form face-to-face contacts. 

The requirement that the particles have the same orientation is globally consistent with a 
lattice packing. Indeed, in the optimal lattice packings of the centrally symmetric Platonic 
solids, each particle has the maximum number of face-to-face contacts that could be obtained 
without violating the impenetrability condition. It is highly unlikely that such particles 
possessing three equivalent principle directions and aligned in the same direction could form 
a more complicated non-lattice periodic packings with densities that are larger than the 
optimal lattice packings. Such non-lattice packings would arise to take advantage of the 
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(a) (b) (c) (d) 

FIG. 13: (color online). Illustration of our conjecture concerning the optimal packings of centrally- 
symmetric particles using two-dimensional octagonal packings. It is shown how different particle 
orientations and arrangements affect the number of face-to-face contacts and the area of the convex 
hulls joining the centers of the neighboring particles around a central one. (a) The octagons do 
not have the same orientation and the number of face-to-face contacts is small, (b) The octagons 
are aligned with 4 contacts per particle to form a lattice packing, but the area of the convex 
hull associated with the central particle is not minimized, (c) The octagons are aligned with 6 
contacts per particle to form a lattice packing, but although the area of the convex hull associated 
with the central particle is smaller than in (b), it is not minimized, (d) The octagons are aligned 
with 6 contacts per particle to form a lattice packing, and the area of the convex hull associated 
with the central particle is minimized. This minimization corresponds to finding the the minimal 
circumscribing hexagon and therefore corresponds to the optimal lattice packing 4l|. 



rotational degrees of freedom. By requiring particle alignments, only translational degrees 
of freedom remain and hence optimization over these degrees of freedom would lead to the 
globally optimal packings, which should be lattice packings. Constraining rotational degrees 
of freedom in this way in a non-lattice packing would at best lead to a local optimum in 
density. This conclusion is clearly supported by our simulations, which use multiple-particle 
configurations in the fundamental cell and only produce the optimal lattice packings. 

In the aforementioned arguments, a key step that is difficult to prove is the observation 
that alignment of each of the particles maximizes the number of face-to-face contacts and 
thus the fraction of space covered by the particles within the convex hull joining the centroids 
of the contacting neighbors around a central polyhedron. A rigorous verification of this step 
would lead to the consistency of the aforementioned local and global optimization criteria. 
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It is noteworthy that in two dimensions, our local criterion to achieve the optimal packings 
of centrally symmetric regular polygons (the two-dimensional counterparts of the centrally 
symmetric Platonic solids) amounts to the identification of optimal neighbor configurations 
with the maximal edge-to-edge contacts per particle such that the convex hull joining the 
centroids of the contacting neighbors is minimal. This statement is rigorously true in two 
dimensions because it is tantamount to a theorem due to Fejes Toth 4l|, which states that 
the densest packing of a centrally symmetric particle in two dimensions can be obtained by 
circumscribing the particle with the minimal centrally symmetric hexagon, which of course 
tiles M^. Finding the minimal circumscribing centrally symmetric hexagon is equivalent to 
finding the minimal area convex hull joining the centroids of the six contacting particles 
(the maximal number can be obtained in two dimensions for centrally symmetric particles). 
The fact that local optimality is consistent with the global optimality in two dimensions 
(e.g., a centrally symmetric hexagon can always tessellate the space) does not hold in three 
dimensions. Although we have presented a three-dimensional generalization of the Fejes 
Toth theorem, by replacing the minimal hexagon with the minimal convex-hull volume, as 
pointed out before, it is extremely difficult to provide a rigorous proof for it. 

Although Conjecture 1 applies to the centrally symmetric Platonic solids, all of our argu- 
ments apply as well to each of the centrally symmetric Archimedean solids. It cannot be true 
for the truncated tetrahedron, which is the only non-centrally symmetric Archimedean solid. 
Specifically, it immediately follows from the results of Ref. Il3| that a non-lattice packing of 
truncated tetrahedra can be constructed based on the "primitive Welsh" tessellation (i.e., by 
removing the small regular tetrahedra) that possesses the density cf) = 23/24 = 0.958333 . . ., 
which is appreciably larger than the optimal-lattice-packing density of 0max = 0.680921 . . ., 
and contains two particle centroids per fundamental cell. In fact, given that the truncated 
tetrahedra cannot tile space, the density of the "primitive Welsh" packing of truncated tetra- 
hedra is so large that it may be the optimal packing of such particles. The lattice vectors 
of this periodic packing are given in Appendix C. 

Since the arguments used to justify Conjecture 1 apply equally well to the 12 centrally 
symmetric Archimedean solids, we are led to the following more general conjecture: 
Conjecture 2: The densest packings of the centrally symmetric Platonic and Archimedean 
solids are given by their corresponding optimal lattice packings. 

The aforementioned arguments can also be extended to the case of superballs, but here 
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the local principle curvatures at the contacting points should be sufficiently small so as to 
maximize the fraction of space covered by the particles within the convex hull joining the 
centroids of the neighbors. The central symmetry and equivalence of the three principle 
directions of superballs means that dense packings of such objects are favored when the 
particles are aligned, which again leads to consistency between local and global optimal- 
ity. The fact that the optimal lattice packing densities of superballs is relatively close to 
the upper-bound values, at least around the sphere point, as well as results from previous 
molecular dynamics simulations also strongly suggest that the densest packings of these 
particles are given by their corresponding optimal lattice packings. 

It is noteworthy that under the assumption that Conjecture 2 is valid, one has upper 
bounds that are the complete analog of f[T^ . i.e., 

(25) 

where Vc is the volume of an appropriately chosen centrally symmetric Platonic or 
Archimedean solid and (p^ax is the corresponding optimal density for such a solid, which ac- 
cording to Conjecture 2 is the optimal lattice packing. This bound will generally be sharper 
than bound f|T^ because the reference optimal packing is less symmetric than the sphere. 
For example, the conjectured bound will be sharp for slightly deformed Platonic and 
Archimedean solids or nonspherical particles derived by smoothing the vertices, edges and 
faces of polyhedra. 

Our work also naturally leads to another conjecture: 
Conjecture 3: The optimal packing of any convex, congruent polyhedron without central 
symmetry generally is not a (Bravais ) lattice packing. 

In other words, the set of such polyhedra whose optimal packing is not a lattice is overwhelm- 
ingly larger than the set whose optimal packing is a lattice. We have seen that because the 
regular tetrahedron and truncated tetrahedron lack central symmetry, dense packings of 
such objects favor face-to-face contacts. Such orientations immediately eliminate the possi- 
bility that lattice packings (in which particles must have the same orientations) are optimal. 
Similarly, it is very plausible that dense packings of most convex, congruent polyhedra with- 
out central symmetry are facilitated by face-to-face contacts and hence the optimal packings 
cannot be lattices. For example, consider a square with one missing corner, i.e., an isosceles 
triangle with a right angle (see Fig. [T^ . At first glance, one might surmise that if the missing 
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(a) (b) 

FIG. 14: (color online). Portions of two packing configurations of pentagons obtained by cutting 
off a corner (isosceles triangle with a right angle) of a square, (a) The optimal lattice packing and 
(b) a two-particle basis periodic packing that tiles the plane. Let side length of the square be 1 
and the lengths of the equal sides of the isosceles triangle he (3 £ (0, 1). The lattice vectors of the 
optimal lattice packing are = i — ^j, eg = (1 — f)i + (l — /3)j and the lattice vectors of the 
periodic packing are ef = i + e^^ = (1 — /3)i + (2 — /3)j with one particle at origin and the other 
at bi = (1 — /9)i + (1 — /3)j, where i, j are the unit vectors along the two orthogonal coordinate 
directions, coinciding with two orthogonal sides of the square. The density (covering fraction) of 
the optimal lattice packing is (l)max = the density of the periodic packing is 

't'max = 1- It can be seen that for all < /? < 1, (j)^ax is always smaller than (f^max- 



piece is sufficiently small, the original lattice packing should still be optimal or nearly opti- 
mal (see Fig. [T^). since lattice packings are optimal for squares. However, no matter how 
small the missing piece may be, a periodic packing in which the fundamental cell contains 
two pentagons can be constructed that tile the plane (see Fig. [Hb). This is done by taking 
advantage of the asymmetry of the particle. Thus, we see from this counterintuitive example 
that if the particle does not possess central symmetry, it is possible to exploit its rotational 
degrees of freedom to yield a periodic packing with a complex basis that are generally denser 
than the optimal lattice packing. On the other hand, there are special cases where the lat- 
tice will be optimal, such as for the rhombic dodecahedron that has one corner clipped j^sj . 
However, these special cases are overwhelmed in number by those whose optimal packings 
are not lattices. If Conjecture 3 is valid, it also appUes to nonspherical particles derived by 
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smoothing the vertices, edges and faces of polyhedra provided that the particle curvature at 
face-to-face contacts is sufficiently small. 

It should not go unnoticed that the densest packings of all of the Platonic and 
Archimedean solids reported here as well as the densest known packings of superballs 
\v\ and ellipsoids 12 1 have densities that exceed the optimal sphere packing density 
4'max = vr/vTs = 0.7408. . .. These results are consistent with a conjecture of Ulam who 
proposed, without any justification, that the optimal density for packing congruent spheres 
is smaller than that for any other convex body 42|. The sphere is perfectly isotropic with 
an asphericity 7 of unity, and therefore, as noted earlier, its rotational degrees of freedom 
are irrelevant in affecting its packing characteristics. On the other hand, each of the afore- 
mentioned convex non-spherical particles break the continuous rotational symmetry of the 
sphere and thus its broken symmetry can be exploited to yield the densest possible packings. 
However, broken rotational symmetry in and of itself may not be sufficient to satisfy Ulam's 
conjecture if the convex particle has a little or no symmetry. 

It will also be interesting to determine whether our conjecture can be extended to other 
polyhedral packings. The infinite families of prisms and antiprisms 4^ provide such a class 
of packings. A prism is a polyhedron having bases that are parallel, congruent polygons and 
sides that are parallelograms. An antiprism is a polyhedron having bases that are parallel, 
congruent polygons and sides that are alternating bands of triangles. Prisms with an even 
number of sides and antiprisms are centrally symmetric. Although prisms and antiprisms 
are naturally grouped with the Archimedean solids (i.e., they are polyhedra in which the 
same regular polygons appear at each vertex), they are generally much less symmetric than 
either the Platonic or Archimedean solids. Moreover, even the centrally symmetric prisms 
and antiprisms generally do not possess three equivalent directions. Thus, it is less obvious 
whether Bravais lattices would still provide the optimal packings for these solids, except for 
prisms that tile space (e.g., hexagonal prism or rhombical prisms). In future work, it would 
be desirable to test whether our conjecture extends to prisms and antiprisms that possess 
central symmetry and three equivalent directions using the ASC scheme. 

It is worth noting that in four dimensions, the analogs of the tetrahedron, cube, octahe- 
dron, dodecahedron and icosahedron are the four-dimensional regular simplex, hypercube, 
cross polytope, 120-cell and 600-cell, respectively. All of these four-dimensional polytopes 
possess central symmetry, except for the simplex 45|]. While the hypercube and cross poly- 
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tope tile W^, the optimal packings of simplices are still likely to be non-lattices. Since our 
conjecture for the three-dimensional Platonic solids should still apply in four dimensions, 
the densest packings of the 120-cell and 600-cell could be their corresponding optimal lattice 
packings. The cross polytopes for d > 5 no longer tile space, and their optimal packings 
may still be their densest lattice packings provided that d is sufficiently small. However, in 
sufficiently high dimensions, the densest lattice packings of the centrally symmetric poly- 
topes are probably no longer optimal, since lattice packings in high dimensions are known 
to possess huge "holes" into which additional particles can be inserted, yielding higher pack- 
ing densities with possibly non-lattice arrangements. Indeed, it has recently been argued 
that disordered sphere packings in very high dimensions could be denser than any ordered 



packing 



46|. 



le production of a wide 



48( 1 ■ icosahedra [49| and 



Recent progress in particle synthesis methods have enabled t 
spectrum of nanoparticle shapes such as tetrahedra |47|, cubes 

prisms js^. In such applications, it would be of great interest to predict the corresponding 
crystal structures, which might possess unusual symmetries and properties. The idea of 
incorporating collective particle motions due to the adaptive fundamental cell should still 
make it efficient to search the desired crystal structures formed by those polyhedral nano- 
building blocks. Finally, it is worth noting that the ASC algorithm is also suitable to generate 
random packings of polyhedral particles. Crucial dynamical parameters of the system, such 
as the strain rate can be properly controlled to produce pac kings with varying degrees of 
disorder, including the maximally random jammed ones 5l||- We will explore disordered 
packings in future work. 
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APPENDIX A: OPTIMAL LATTICE PACKINGS OF THE PLATONIC AND 
ARCHIMEDEAN SOLIDS 



Here we collect fundamental packing characteristics (i.e., the lattice vectors, densities and 



contact numbers) of the optimal lattice 
most of which are obtained from Refs. 



mc 



22 



dngs of the Platonic and Archimedean solids. 



34 and 



23l . The Platonic solids are explicitly 



defined as the regions (sets of points) bounded by a set of linear equations of the coordinates. 
The Archimedean solids are either defined as intersections of different Platonic solids or 
delineated based on their symmetry. A point set 5* representing the polyhedron multiplied 
by a number a means an isotropic expansion of the polyhedron 5" with linear ratio a. n5'2 
represents the intersection region of two polyhedra S\ and 5*2. The lattice vectors are given 
by column vectors (i.e., the basis vectors are the unit vectors along the coordinate axis). 
The naming code used below for the Platonic and Archimedean solids is the same one used 
in Figs. 1 and 3. 

A tetrahedron (PI) has 4 vertices, 6 edges and 4 triangular faces. It is defined as the 
region 

PI = {x G : Xi + a;2 + < 1, —X\ — 2:2 + < 1, —X\ + 2:2 — 0:3 < 1, Xi — X2 — 0:3 < ij . 

(Al) 

The optimal lattice vectors are given by 

ai = (2, -1/3, -1/3)'^ , a2 = (-1/3, 2, -1/3)'^ , = (-1/3, -1/3, 2f , (A2) 

where the superscript "T" denotes transpose of a column vector. Each tetrahedron of the 
packing contacts 14 others. The packing density is 4>^ax — 18/49 = 0.367346 . . .. 

An icosahedron (P2) has 12 vertices, 30 edges and 20 triangular faces. It is defined as 
the region 

P2 = {x eM? : \xi\ + |X2| + |X3| < 1, |<l>Xi| + |X3/$| < 1, |<I>X2| + |xi/<I>| < 1, |<I>X3| + |X2/$| < 

(A3) 
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where $ 



ai(x) 



'1 + is the golden ratio. The lattice vectors are given by 

( (-33/8 - 39y5/8)x2 + (39/4 + 33V5/4)x - 11/4 - 3v^/2 ^ 
-1/4 - y/h/A)x + 1 + v^/2 

-19/2 - 8V5)x + 13/4 + 3^5/2 



(1+*) 



(33/8 + 39\/5/8)x2 



ai2[x) 



(i+i>) 



( (-39/8 - 33v/5/40)x2 + (35/4 + 41v^/20)x - 5/2 - 23v^/20 ^ 
(5/4 + \/h/A)x - 1 - V5/2 
(-39/8 - 33V5/40)x2 + (15/2 + 9v^/5)x + 13/4 - 3v^/20 



(A4) 



( (3/2 + v^/2)x - 2 - V5 ^ 



as ((a;)) 



(1+*) 



v 



X 





where x G (1, 2) is the unique root of the polynomial 

1086x^ - (1063 + 113v/5)a;2 + (15^5 + 43)x + 102 + 44^5 = 0. 



(A5) 



It is found that x = 1.59160301 . . . and therefore the lattice vectors, up to nine significant 
figures, are given by 



ai = (0.711782425,0.830400102,1.07585146)' 
a2 = (-0.871627249,0.761202911,0.985203828) 



(A6) 



ag = (-0.06919791,1.59160301,0) 



T 



max 



Each icosahedron of the packing contacts 12 others. The packing density is 
0.836357... 

A dodecahedron (P3) has 20 vertices, 30 edges and 12 pentagonal faces. It is defined as 
the region 



P3 = {a; G : + jxal < 1, I^Xal + |xs| < 1, |$Xs| + < l} 



The optimal lattice vectors are given by 



(A7) 



ai = (2/(1 + $), 2/(1 + <l>),0)' 
a2 = (2/(1 + $), 0,2/(1 + $)) 



T 



(A8) 



as = (0,2/(1 + $), 2/(1 + $))^. 
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Each dodecahedron of the packing contacts 12 others. The packing density is 4>^ax = 
(2 + $)/4 = 0.904508.... 

An octahedron (P4) has 6 vertices, 12 edges and 8 triangular faces. It is defined as the 
region 

P4 = {x e : |xi| + Ixal + \xs\ < l} . (A9) 
The optimal lattice vectors are given by 

ai = (2/3, 1, 1/3)'^ , a2 = (-1/3, -2/3, 1)'^ , ag = (-1, 1/3, -2/3)'^ . (MO) 

Each octahedron of the packing contacts 14 others. The packing density is (p^ax — 18/19 = 
0.947368.... 

A cube (P5) has 8 vertices, 12 edges and 6 square faces. It is defined as the region 

P5 = {x e : < 1} . (All) 

The optimal lattice vectors are given by 

ai = (2,0,0)^, a2 = (0,2,0)^, as = (0,0, 2)^. (A12) 

Each cube of the packing contacts 26 others (which includes vertex-to- vertex contacts). The 
packing density is (p^ax = 1- 

A truncated tetrahedron (Al) has 12 vertices, 18 edges and 8 face: 4 hexagons and 4 
triangles. It is defined as the region 

Al = {x eR^ : X e 5- Pin-3 ■ Pi} . (A13) 

The optimal lattice vectors are given by 

ai = (4/3, 4, 8/3)^ , a2 = (4, -8/3, -4/3)^ , ag = (-8/3, 4/3, -4f . (A14) 

Each truncated tetrahedron of the packing contacts 14 others. The packing density is 
= 207/304 = 0.680921 . . .. 
A truncated icosahedron (A2) has 60 vertices, 90 edges and 32 faces: 20 hexagons and 12 
pentagons. It is defined as the region 

A2 = {x G : X G (1 + $) ■ P2 n (4/3 + $) ■ P3} . (A15) 
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The optimal lattice vectors are given by 

ai = (1 + <l>)af2, as = (1 + $)af , ag = (1 + <l>)aP, (A16) 

where af^ (i = 1,2,3) are the lattice vectors of the optimal lattice packing of icosahedra. 
Each truncated icosahedron of the packing contacts 12 others. The packing density is 
= 0.784987... 

A snub cube (A3) has 24 vertices, 60 edges and 38 faces: 6 squares and 32 triangles. Let 
the snub cube orientated in a way such that its 6 square faces lie in the hyperplanes 

{x eR^ ■.\xi\ = l{i= 1,2,3)} . (A17) 

The optimal lattice vectors are given by 

ai = (2,0,0)^, a2 = (0,0,2)^, ag = (1, 2/y* - 2, -1)^ . (AlS) 

where y* is the unique real solution of + + y = 1. Each snub cube of the packing 
contacts 12 others. The packing density is 0^^^ = 0.787699 . . . 

A snub dodecahedron (A4) has 60 vertices, 150 edges and 92 faces: 12 pentagons and 80 
triangles. Let the snub dodecahedron orientated in a way such that its 12 pentagonal faces 
lie in the hyperplanes of the faces of the dodecahedron (1 + $) ■ P3. The optimal lattice 
vectors are given by 

ai = (2,2,0)'^, a2 = (2,0,2)'^, a3 = (0,2,2)^. (A19) 

Each snub dodecahedron of the packing contacts 12 others. The packing density is (p^^x ~ 
0.788640 . . . 

A rhombicosidodecahdron (A5) (also known as small rhombicosidodecahedron) has 60 
vertices, 120 edges and 62 face: 12 pentagons, 30 squares and 20 triangles. It is defined as 
the region: 

v45 = ja; G : a; G (3$ + 2) ■ [P5 n (1 + V2)PA] D (4$ + 1) ■ P2 n (3 + 3$) ■ P3} . 

(A20) 

The optimal lattice vectors are given by 

ai = (($ - 1)/(2<I> + 1), 7, (9$ + 4)/(2<l> + 1))^ , 

a2 = ((9$ + 4)/(2$ + 1), ($ - 1)/(2<I> + 1), 7f , (A21) 

as = (7, (9$ + 4)/(2$ + 1), (<f - 1)/(2<I> + 1))^ . 
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Each rhombicosidodecahedron of the packing contacts 12 others. The packing density is 
0^^^ = (8$ + 46)/(36$ + 15) = 0.804708 . . .. 

A truncated icosidodecahdron (A6) (also known as great rhombicosidodecahedron) has 
120 vertices, 180 edges and 62 faces: 12 decagons, 20 hexagons and 30 squares. It is defined 
as the region 

= |x e : X G (5$ + 4) ■ [P5 n (1 + V2)PA] n (6$ + 3) ■ P2 n (5 + 5$) ■ P3| . 

(A22) 

The optimal lattice vectors are given by 

ai = (10,10,0)^, a2 = (10,0,10)^, ag = (0, 10, 10)^ . (A23) 

Each truncated icosidodecahedron of the packing contacts 12 others. The packing density 
is 0^^^ = (2<l>/5 + 9/50) = 0.827213 . . .. 

A truncated cubeoctahedron (A7) (also known as great rhombicubeoctahedron) has 48 
vertices, 72 edges and 26 faces: 6 octagons, 8 hexagons and 12 squares. It is defined as the 
region 

A7 = {a; e . i^,^] ^ i^,^! < (2 + 3V2), \x2\ + jxal < (2 + 3^2), 
^2! + kal < (2 + 3v^),x G (2v^+ 1) ■ P5 n (3v^ + 3) • P4}. 
The optimal lattice vectors are given by 

, T 
, T 



(A24) 



ai = (4v^+2,-4V2-l + 2a,4V2 + l-2a)^ 
a2 = {V2/2 -3/2 + a, -3v^/2 + 1/2 + a, 4v^ + 2) ' , (A25) 
ag = (7/2 + 7^2/2 - a, 1 + 2a, 5^2/2 + 3/2 - af 

where a = y/33{V2 + l)/6. Each truncated cubeoctahedron of the packing contacts 12 
others. The packing density is ip^ax = 0.8493732 . . .. 

An icosidodecahedron (A8) has 30 vertices, 60 edges and 32 faces: 12 pentagons and 20 
triangles. It is defined as the region 

A8= {x eM.^ : X e P2n P3} . (A26) 

The optimal lattice vectors are given by 

ai = (2/(1 + $), 2/(1 + $),0)^, 

a2 = (2/(1 + $), 0,2/(1 + $))^, (A27) 
a3 = (0,2/(1 + $), 2/(1 + $))^. 
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Each icosidodecahedron of the packing contacts 12 others. The packing density is 4>max = 
(14 + 17<l')/48 = 0.864720.... 

A rhombicuboctahedron (A9) (also known as small rhombicubeoctahedron) has 24 vertices, 
48 edges and 26 faces: 18 squares and 8 triangles. It is defined as the region 

y49 = {x e : Ixil + |x2| < 2, |x2| + bal < 2, 

I II I I _ , I I I Ji - , ^^28) 

kil + lxal <2,xey2-P5n(4- V2)-P4}. 
The optimal lattice vectors are given by 

ai = (2,2,0f, a2 = (2,0,2f, a3 = (0,2,2f. (A29) 

Each truncated cubeoctahedron of the packing contacts 12 others. The packing density is 
0^^^ = (16v^ - 20)/3 = 0.875805 . . .. 

A truncated dodecahedron (AlO) has 60 vertices, 90 edges and 32 faces: 12 decagons and 
20 triangles. It is defined as the region 

AW = {x G (1 + $) ■ P3 n [(7 + 12<l>)/(3 + 4$)] ■ P2} . (A30) 

The optimal lattice vectors are given by 

ai = (2,2,0)^, a2 = (2,0,2)^, ag = (0, 2, 2f . (A31) 

Each truncated dodecahedron of the packing contacts 12 others. The packing density is 
<PL. = (5$ + 16)/(24<|. - 12) = 0.897787. . .. 

A cuboctahedron (All) has 12 vertices, 24 edges and 14 faces: 6 squares and 8 triangles. 
It is defined as the region 

All = {x G P5n2-P4}. (A32) 
The optimal lattice vectors are given by 

ai = (2, -1/3, -1/3)^ , a2 = (-1/3, 2, -1/3)^ , ag = (-1/3, -1/3, 2f . (A33) 

Each cuboctahedron of the packing contacts 14 others. The packing density is 4>^ax — 
45/49 = 0.918367.... 

A truncated cube (A12) has 24 vertices, 36 edges and 14 faces: 6 octagons and 8 triangles. 
It is defined as the region 

A12 = {x G P5 n (1 + V2) ■ P4} . (A34) 
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The optimal lattice vectors are given by 

ai = (2,-2«,0f , a2 = (0,2,-2«f , as = (-2«, 0, 2f , (A35) 

where a = (2 — -\/2)/3. Each truncated cube of the packing contacts 14 others. The packing 
density is (p^^^ = 9/(5 + 3^2) = 0.973747 . . .. 

A truncated octahedron (A13) has 24 vertices, 36 edges and 14 faces: 8 hexagons and 6 
squares. It is defined as the region 

^13 = {x e P5 n (3/2) ■ P4} . (A36) 
The optimal lattice vectors are given by 

ai = (2,0,0f, a2 = (2,2,0f, a3 = (l,l,-lf. (A37) 

Each truncated octahedron of the packing contacts 14 others. The packing density is (p^ax = 
1. 



APPENDIX B: LATTICE VECTORS AND OTHER CHARACTERISTICS OF 
THE DENSEST KNOWN TETRAHEDRAL PACKING 

Here we report the lattice vectors of fundamental cell for the densest known tetrahedral 
packing with 314 particles up to 12 significant figures, although our numerical precision 
is not limited by that. The side length of the tetrahedron is do = 1.8 and its volume is 
Vp = V2dl/12 = 0.687307791313. The lattice vectors are given by 

ai = (6.130348985438, - 1.714011427194, 0.022462185673)^, 
as = (1.622187720300, 6.106123418234, 0.123567805838)^, (Bl) 
ag = (0.055014355972, 0.116436714923, 6.526427521512)^. 

The volume of the fundamental cell is Vol(F) = |ai x aa ■ ag] = 262.344828467328. Thus, 
the packing density is readily computed as 

N ■ v„ 314 X 0.687307791313 , , 

^- P - -. 0.822637319490. (B2) 



Vol(F) 262.344828467328 



The coordinates of each of the 314 tetrahedra are given elsewhere 371]. 
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APPENDIX C: LATTICE VECTORS OF THE DENSEST KNOWN PACKING 
OF TRUNCATED TETRAHEDRA 



The "primitive Welsh" tessellation of space consists o 



dra and small regular tetrahedra, as described in Ref. 



13 



" truncated large regular tetrahe- 
When the small tetrahedra are 



removed, the remaining truncated tetrahedra give the densest known periodic (non-lattice) 



packing of such objects. Following Ref. 



13, the centroids of the truncated tetrahedra sit 



at the nodes and 1 of the two adjacent cells of a body-centered cubic lattice shown in 
Fig. 2(a) of Ref. [l3|. Let length of the edges of the truncated tetrahedron be (Ie = "\/2/2, 



then the lattice vectors are given by 



ai = (1, 1, 0)'^, a2 = (1, 0, 1)^, ag = (0, 1, 1)'^. 



(CI) 



The basis vectors for the centroids of 0-type and 1-type truncated tetrahedra are respectively 
given by 

bo =(1/2, 1/2, 1/2)^, bi = (0,0,0)^. (C2) 
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